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Abstract 

We analyze the effect of vr^-terms in the QCD perturbative expansions for the 

(^ ' s-channel effective coupling and observables, the effect known from the 80s. We 

^^ . remind that these terms can be collected into specific functions — strong s-channel 

coupling d(s) and its effective powers 2lfc(s) free of ghost singularities. Further on, 

we study the structure of perturbation theory for observables and its reformulation 

in terms of nonpower perturbation expansion over the set {2lfc(s)} . 

cn ! Then we discuss the influence of this effect on the numerical values of a^ as 

extracted from experiments. The main result is that the common two- loop (NLO, 

^ ■ NLLA) approximation widely used in the five-quark (10 GeV < ^/s < 170 GeV) 

K*" , region for a shape analysis contains a systematic negative error of a 1-2 per cent 

^^ \ order of magnitude for the extracted cks ■ Our physical conclusion is that the 

as{M'^) value averaged over the / = 5 data 

C^ ! 

g : < as{M^) >f=5^ 0.124 



(N 



appreciably differs from the currently accepted "world average" (= 0.118). 



o 

(D : 1 Preamble 

^^ 

Usually, physical quantities in the time-like channel, like the cross-section 
K> ! ratio of the inclusive e^e^ -^ hadron annihilation or the r-decay process, 

H I are presented in the form of two- or three-term perturbation expansion 



Rq 



l + r{s); r{s) = cias{s) + C2ag + C3ag + . . . (1) 



(our coefficients Ck = Ck vr~^ are normalized differently from the commonly 
adopted, like in Refs.|l|, 0, |^) over powers of effective QCD coupling a^ which 
is supposed ad hoc to be of the same form as in the Euclidean domain, e.g.. 



^""^" - ^-l^ + ^I'^c-'^-'-^-i'+^^l^ 



a 



+ 



1 



h\ ( -\yi^ L + -\yi^ L + 2hiL --\ -36162 lnL + - 



Here, L = ln(s/A^) and for the beta-function we use normalization 

/3(a) = -/5o a^ - /9i a^ - /Sa a"^ + . . . = -/So a^ (l + 61 a + 62 a^ + . . .) 
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that is also free of it powers. Numerically, 

33-2/ ^ ,^, 153-19/ 



W) 



12tt 



biif) 



; &i(4±l) = 0.490: 



-0.089 
h0.076 



Coefficients Ck>3 
lower Ck'- 

, (^/?o(/))'ci 



27r(33-2/) 
dk — Sk include "tt^ structures" 6k proportional to 



<54 = (^/?o)Mc2 + -6iCi); vr^/3o'(4±l) 



4M0-fi . 



(2) 



These structures 6k arise 0, ^, ||, ^ in the course of analytic continuation 
from the Euclidean to Minkowskian region. Coefficients d^ should be treated 
as a genuine fcth-order ones. Just they have to be calculated with the help 
of relevant Feynman diagrams. 

To illustrate, consider the three-flavor case for r-decay, / = 4 , 5 cases 
for e+e^ —>■ hadron annihilation and Zq decay (with / = 5) — see Table 1 in 
which we also give values for the vr^-terms . 

Table 1 



Process 


f 


Cl 


C2 = ^2 


C3 


ds = C3- 63 


^3 


^4 


r decay 


3 


l/vr 


.526 


0.852 


1.389 


0.537 


5.01 


e"'"e~ 


4 


.318 


.155 


-0.351 


0.111 


0.462 


2.451 


e"'"e~ 


5 


.318 


.143 


-0.413 


-0.023 


0.390 


1.752 


Zq decay 


5 


.318 


.095 


-0.483 


-0.094 


0.390 


1.576 



Here, all coefficients Ck , dk and 6k , due to normalization ([|), are of an 
order of unity. One can see that, in the high energy region, contribution of 
63 prevails in C3 . 



2 Preliminary quantitative estimate 

In practice, the vr^-terms often dominate in higher expansion coefficients. 
This effect is especially strong in the / = 5 region. Meanwhile, just in 
this region people often use the so-called NLLA approximation, that is the 
two-term representation 



0(s)=Ci(a,/7r) + C2(«,/7r) 



(3) 



for an observable 0{s) when next, the three-loop, coefficient C3 is not known. 
This is the case, e.g., with event-shape § analysis. 

On the basis of the numerical estimates of Table 1, in such a case, we 
recommend to use the three-term expression 



0^{s) = rfi I 



a^ — 



^'PLs 



aA + d2a. = ciUs + C2a. - 6s a^ 



(4) 



i.e., to take into account the known predominant vr^ part of the next coeffi- 
cient C3 . As it follows from the comparison of the last expression with the 
previous, two-term one, the a^ numerical value extracted from eq.(^), for 
the same measured value Oobs, will differ by a positive quantity (e.g., in the 
/ = 5 region with a^ ~ 0.12 -^ 0.15) 






1 + 2nd2as 



^=^ 1.225 a^ 

= -^^ - 0.002 - 0.003 

20-MooGcV ^ + ^.90 as 



that turns to be numerically important. 

Moreover, in the / = 4 region, where the three-loop approximation is 
commonly used in the data analysis, the vr^ term ^4 of the next order turns 
out also to be essential. Hence, we propose to use the four-term expression 



O4 (s) = dias + d2 a^ + C3 a^ - 64 a^ ; C3 = d^ - 63 (5) 

(instead of the three-term one (|I|)) that is equivalent to 

of (s) =dila,- -^a^ -h- T^^l a\Wd2 [o?, - t^^I a\] + da o?, (6) 

with (^3 and b^ defined^, in eq.(|^). 

The three- and two-term structures in curly brackets are related to spe- 
cific expansion functions a and 21 defined below (|TD|) and entering into the 
non-power expansion (|TTD . 

To estimate roughly the numerical effect of using this last modified expres- 
sion (^, we take the case of e"'"e~ inclusive annihilation. For y^ ~ 3-i-5GeV 
with a, ^ 0.28 ^ 0.22 one has 






1 + lud^as 



^=^ 1.07 a^ 

= ' ~ 0.005 - 0.002 

3^5GeV 1 + 0.974 a, 



— an important effect on the level of ca 1 -h 2% . 

Moreover, the (Aas)4 correction turns out to be noticeable even in the 
lower part of the / = 5 region! Indeed, at ^/s ~ 10 -i- 40 GeV with a, ~ 
0.20 ^ 0.15 we have 

(Aa,)4|fot40GeV ^ 0.71 a\ ^ (1.1 - 0.3) ■ 10"^ (< 0.5%) . 

3 Non-power expansion in the Minkowskian region 

The so-called vr^ terms in the s-channel perturbative expansions for the 
invariant coupling and observables have a simple origin. 



As it is well known, the usual invariant coupling originally defined 
in terms of real constants Zi, counter-terms of finite Dyson renormalization 
transformation, can be expressed via a product of dressed symmetric vertex 
and propagator amplitudes taken at space-like values of their arguments. 



a{Q^a)=ar\Q\a)Y[di{Q' 



a) 



Hence, by construction, it is a real function defined in the Euclidean region. 
Transition to the time-like region, with logs branching InQ"^ -^ Ins — m 
transforms all relevant amplitudes into complex functions r(s, a), (ij(s, a) . 
Here, the problem of appropriate defining of effective coupling in the time- 
like domain arises. 

For this goal, we shall follow the idea devised in the early 80s by Radyushkin 
1^ and Krasnikov-Pivovarov |]^ . There, an integral transformation R reverse 
to the dipole representation for the Adler function has been used. 

We propose to treat this representation as an integral operation 



R{s) ^ D{z) = g^ 



ds 



is + z) 



R{s) =D{i?(s)} 



(7) 



transforming a function R{s) of a real positive (time-like) argument into a 
function D{z) given in the cut complex plane with analytic properties equiv- 
alent to those following from the Kallen-Lehmann integral representation. In 
particular, the function D{Q'^) is real on the positive (space-like) real axis 
at 2 = Q2 + ^0 ; Q2 > . 

The reverse operation is expressible in the form of a contour integral 



R{s) = - / -D,,{~z)^K[D{Q' 



With the help of the latter, one can define ITT], |I^ an effective invariant 
time-like coupling a{s) = R[q;s(Q^)] . Omitting some technical details, we 
give a few resulting Q, ^, [T^ expressions. 



r,(l) 



-1 



E.g., starting with one-loop al = [/3oln((5 /A )] one has R 



aP 



a^^Hs) = — 



Po 



1 1 L 
arctan — 

2 TT n 



L>0 



1 TT S 

— — arctan — ; L = In — - . (8) 

Poll L A^ 



At the same time, to iai {Q^)) and iai (Q^)) there correspond 



^?is) 



R 



(«i^^) 
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and ^'i\s) 



/3l [L2 + n- 



n2 • 



In the two-loop case, for a "popular" expression 

1 . ,.Anl 



Q' 



J-W)-l^-, ^ = ln^ 



one obtains PI the two- loop "pop" effective s-channel coupling 



'^pop\^) 



1 + 



biL 



L2 



vr" 



aW| 



^)- 



bi In [y/UT^ + 1 



Po 



L^ + tt"^ 



(9) 



Both the expressions (^ and (^ are monotonically decreasing with a finite 
IR (5(0) = l//5o(/ = 3) ^ 1.4 value. Meanwhile, higher functions go to the 
zero 2lfc(0) = at the IR limit. 

In the case L :» vr, it is possible to expand a and 2tfc in powers of tt^/L^. 
Then functions a and 2I2 can be presented as expansions in powers of common 
(is — 1/L. They correspond to curly brackets in (^. 

In [|, Q, as a starting point for observables in the Euclidean, i.e., space- 
like domain Q^ > 0, the perturbation series 

k>l 

has been assumed. It contains powers of usual, RG summed, invariant cou- 
pling asiQ"^) that obeys unphysical singularities in the infrared (IR) region 
around Q^ ^ A| . 

By using the R transformation, we obtain in the Minkowskian region the 
"transformed" expansion over a non-power set of functions 

R^{s) = R [Dpt(g')] = 1 + J2dkMs) ; Ms) = R K(Q')] (10) 

fc>l 

free of the mentioned singularities. Properties of these functions have been 
analyzed in detail in our previous paper [O] — see also Ref. Jlli. For a more 



detailed numerical information on the functions a , 2I2 and 2I3 see Ref. 

Here, we give condensed information that will be enough for a few illus- 
trations. 

Table 2 

Three-loop APT results for A^^ = 290 GeV; ^^(Mf) = 0.125 



Vi/GeV 


5 


10 


15 


20 


30 


50 


60 


90 


150 


as{s) 


.235 


.195 


.177 


.165 


.153 


.137 


.133 


.125 


.115 


d{s) 
IO2I2 
IOO2I3 


.221 
.456 
.871 


.186 
.330 

.555 


.170 
.275 
.436 


.160 
.246 
.357 


.148 
.214 
.299 


.136 
.180 
.232 


.132 
.169 
.213 


.123 
.149 
.177 


.114 
.129 
.143 



Both in the Figure 1 and in Table 2, we give 3- loop solutions for otg as 
well as for the modified, so-called global (for detail, see paper |]l3l) functions 
a = 2ti , 2I2 and 2I3 calculated within the MS scheme for the cases A(5) = 



215 GeV, a,(M: 



0.118 and A(5) = 290 GeV, as(M|) = 0.125. 



0,4-1 



0,3- 



0,2 



0,1 - 




Figure 1: Effective global Minkowskian, a , and Euclidean, 
functions, as compared with the standard one ag fat A (5) = 

o.iis;. 



Oan expansion 
-- 350 MeV and 



o^siMl) 



We have chosen these two cases as limiting ones as far as in many practical 
cases real figures lie between these limits. 

In the first figure we give three curves a^ , a and a^n related to the same 
physical case for A3 = 350 MeV and as^M"^) = 0.118. The curves a and ttan 
on the figure go a bit slanting than usual, the ag , dotted curve. This is quite 
natural, as they both are regular in the vicinity of the A singularity. 

Meanwhile, only two first, a and Oan have direct physical meaning (com- 



pare with conclusion of |T3[)- Just their values have to be determined from 
any given experiment. Nevertheless, in the four- and five-flavour regions one 
can still refer to a^ and as{M'^) as to traditional theoretical objects. 
Now, instead of (|l]), with due account to ([10|) , we have 



r{s) 



a[s\ 



TT 



+ ^2212(5) + 4 213(5) 



(11) 



with beautifully decreasing coefficients dk . Just this nonpower expansion, 
strictly speaking, should be used instead of its approximations, eqs.(|D and 
(^), for data analysis in the time-like region. 

At the same time, in the Euclidean, we have also non-power expansion 

d(Q^) = ^d9!l + d, A(Q') + ds AsiQ') (12) 

IT 

that can be related to (|ll|) by transformation (|^ in the framework of Invari- 
ant Analytic Approach (refs.|jT6|, p!7[|). 

These non-power expansions, free of unphysical singularities, jointly form 
a correlated system. The latter has been studied in detail in Refs. ||13[ and 



T8| . We call it Analytic Perturbation Theory (APT). 



4 Numerical illustrations 

To illustrate, let us start with a few cases in the / = 5 region. 

To begin with, consider the T decay. According to the Particle Data 
Group (PDG) overview (see their Fig.9.1 on page 88 of Ref. |]l|), this is (with 
as(Mx) — 0.170 and as(M|) = 0.114) one of the most "annoying" points 
of their summary of as(M|) values. It is also singled out theoretically. The 
expression for the ratio of decay widths starts with the cubic term 

R{T) = Roal{Mr)il + eias) with d ^ 1 . (13) 

Due to this, the vr^ correction^ is rather big here 

2t3 ^ «^ (1 - 2{7Tf3ofa!) . (14) 

Accordingly, 

Aas{Mr) = ^ {nPo)^ al{Mr) ^ 0.0123, 
that corresponds to 

Aa^(M^) = 0.006 with a,(M^) = 0.120 . (15) 

Now, let us turn to a few cases analyzed by the three-term expansion 
formula (|^). For the first example, take e~^e~ hadron annihilation at i/i = 
42GeVandllGeV. 

A common form (see, e.g., Eq.(15) in Ref.p|) of theoretical presenting of 
the QCD correction in our normalization looks like 

r,+,-(s) = 0.318a,(s) + 0.143a^- 0.413 a^ (16) 



^First proposal of taking into account this effect in the T decay was discussedg more 
than a quarter of century ago. Nevertheless, in current practice it is neglected. 



Starting with re+e-(42) ~ 0.0476, one has as(42) = 0.144. Along with our 
new philosophy, one should use instead 

r^+^~{s) = 0.318 a(s) + 0.143 2t2(s) - 0.023 2t3(s) (17) 

that yields a(42) = 0.142 with as(42) = 0.145 and as{M^) = 0.127 to be 
compared with otsiM"^) = 0.126 under a usual analysis. 

Quite analogously, for re+e-(ll) — 0.0661; as(ll) = 0.200, we obtain 
a(10) = 0.190 that corresponds to asiM'^) = 0.129 instead of 0.130. 

For the next example, we take the Zq inclusive decay. Experimental ra- 
tio Rz = r(Zo -^ hadrons)/T{Zo -^ leptons) = 20.783 ± .029 is usually 
presented as follows: Rz = Ro{l + rziM^)) with Rq = 19.93 . A common 
form (see, e.g., Eq.(15) in Ref.0]) of presenting of the QCD correction in our 
normalization looks like 

rz{Ml) = 0.3326a, + 0.0952 a^ _ 0.483 a^ . 

To [rz]^f,^ = 0.04184 there corresponds a,{Ml) = 0.1241 with A^^ = 
292 MeV . In the APT case,from 

rziMl) = 0.3326 a(M|) + 0.0952 2t2(M|) -0.094 2l3(M|) (18) 

we obtain a(M|) = 0.122 and a,(M|) = 0.124 that relates to A^^) = 
290 MeV . Note that here the three-term approximation of (|^) gives the same 
relation between the as(M|) and a(M|) values. 

Nevertheless, in accordance with our preliminary estimate for the {Aas)4 
role, even the so-called NNLO theory needs some tt^ correction in the W = 
y/s < 50 GeV region. 

Now, turn to the experiments in the HE Minkowskian (mainly with a 
shape analysis) that usually are confronted with two-term expression (|^). As 
it has been shown below, the main theoretical error in the / = 5 region can 
be expressed in the form 

(Aa,(s)|^o=;iooGcV ^ 1.225 a^(s) ^ 0.002 - 0.003 . (19) 

An adequate expression for the shift of an equivalent as{M'^) value is 

[Aa,(M|)]3 = 1.225a,(s)a,(M|)2 . (20) 



We give results of our approximate APT calculations, mainly by Eqs. ([19D 
and (|20|) , in the form of Table 3 and Figure 2. At the last column of the Table 
3 in brackets we indicate difference between the APT and usual analysis. By 
bold figures the results of the three-loop analysis are singled out. 
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Table 3 

The APT revisedQ part (/ = 5) of Bethke's[j2| Table 6 





v^ 


loops 


as (s) 


a,(m^) 


tts (s) 


asim'i) 


Process 


GeV 


No 


ref.[2] 


ref.[2] 


APT 


APT 


T-decay [| 


9.5 


2 


.170 


.114 


.182 


.120 (+6) 


e^e~[ahad] 


10.5 


3 


.200 


.130 


.198 


.129(-1) 


e'^e~[j k.sh] 


22.0 


2 


.161 


.124 


.166 


.127(+3) 


e'^e~[j k.sh] 


35.0 


2 


.145 


.123 


.149 


.126(+3) 


e^e'[ahad] 


42.4 


3 


.144 


.126 


.145 


.127(+1) 


e+e~[j &s/i] 


44.0 


2 


.139 


.123 


.142 


.126(+3) 


e+e~[j &s/i] 


58 


2 


.132 


.123 


.135 


.125(+2) 


Zq — ^ had. 


91.2 


3 


.124 


.124 


.124 


.124 (0) 


e~^e^[j Lsh] 


91.2 


2 


.121 


.121 


.123 


.123(+2) 


e~^e~[j Szsh] 


133 


2 


.113 


.120 


.115 


.122(+2) 


e~^e~[j &s/i] 


161 


2 


.109 


.118 


.111 


.120(+2) 


e~^e~[j &ish] 


172 


2 


.104 


.114 


.105 


.116(+2) 


e~^e~[j ksh] 


183 


2 


.109 


.121 


.111 


.123(+2) 


e~^e~[j &s/i] 


189 


2 


.110 


.123 


.112 


.125(+2) 



Averaged < a^(M2) >j=5 values 0.121; 0.124; 

°"j & sh" — jets and shapes; Figures in brackets in the last column give the 
difference Aas(M^) between common and APT values. 
''Taken from Ref.^l. 



Let us note that our average over events from Table 6 of Bethke's review 
nicely correlates with recent data of the same author (see Summary of 
Ref.jg). The best x^ fit yields as(M|)[2] = 0.1214 and asiM^)APT = 
0.1235 . This gives minimum xL = 0.197 and Xapt ~ 0.144 with impressive 
ratio (~ 0.73) illustrating the effectiveness of the APT procedure. 

On the Fig.2 by open circles and bullets (o, •) we give two- and three- 
loops data mainly from Fig. 10 of paper [^. The only exclusion is the T decay 
taken from the Table 6 of the same paper. By crosses we marked the new 
'APT values" calculated approximately mainly with help of Eq. ( |T9|) . 

For clearness of the vr^ effect, we skipped the error bars. They are the 
same as in the mentioned Bethke's figure and we used them for calculating 



aJs) 



0.18 



0.16 



0.14 



0.12 



A^'s/MeV 



a,(M|) 

0.1235 

0.1184 




Figure 2: The new APT analysis for Og in the five-flavour time-hke region. Crosses 
(+) differ from circles (o,») by vr^ correction (^). Solid APT curve relates to 



A 



(5) 
MS 



270 MeV and 



as{Ml) 



(dot-and-dash curve) Og (at A*^^-* 
Fig. 10 of paper @. 



0.124. To compare, we give also the standard 
= 213 MeV and as{M];) = 0.118) taken from 



5 Conclusion 

We have established a few qualitative effects: 

1. Effective positive shift Aa^ = +0.002 in the upper half (> 50 GeV) of 
the / = 5 region for all time-like events that have been analyzed up to now 
in the NLO mode. 

2. Effective shift Aa, ~ +0.003 in the lower half (10 ^ 50 GeV) of the 
/ = 5 region for all time-like events that have been analyzed in the NLO 
modes. 



3. The new value 



(21) 



as{M'z) = 0.124 

by averaging over the / = 5 region. 

These results are based on a plausible hypothesis on the "tt^- terms" 
prevalence in expansion coefficients for observable in the Minkowskian do- 
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main. The hypothesis has some prehminary support but needs to be checked 
in a more detail. 

Nevertheless, our result ( pT]) being taken as granted, rises two physical 
questions: 

- The issue of self-consistency of QCD invariant coupling behavior be- 
tween the "medium (/ = 3,4)" and "high (/ = 5,6)" regions. 

- The new "enlarged value" (PTf) can influence various physical specula- 



tions in the several hundred GeV region. 
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